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A. J, RICHARDSON 
J. H. EVERITT 
C. L. WIEGAND 
USDA, SEA, AR 

Soil and Water Cons er vat ion Research' 

’ ~ Weslaco, Texas 78596 

Estimating Total Standing Herbaceous Biomass Production With 
LANDSAT MSS Digital Data 

The Landsat derived perpendicular vegetation index may be useful 
to range managers in estimating herbaceous biomass production. 

ABSTRACT: Rangeland biomass data were correlated with spectral vegeta- 

tion indices, derived from LANDSAT MSS data. LANDSAT data from five 
range and three other land use sites in Willacy and Cameron Counties 
were collected on October 17 and December 10, 1975, and on July 31 and 
September 23, 1976. The overall linear correlation of total standing 
herbaceous biomass with the LANDSAT derived perpendicular vegetation 
index was highly significant (r - 0.90**) for these four dates. The 
standard error of estimate was 722 kg/ha. Biomass data were recorded 
for two of these range sites for 8 months (March through October) during 
tka 1976 growing season. Standing green biomass accounted for most of 
the increase in herbage, starting in June and ending about September and 
October. These results indicated that satellite data may be useful for 
the estimation of total standing herbaceous biomass production that could 
aid range managers in assessing range condition and animal carrying 
capacities of large and inaccessible range holdings. 
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INTRODUCTION 

To make proper range management decisions, range managers need to 
know more about (1) the duration of the safe grazing season in relation 
to weather conditions and current herbage production, (2) how to control 
seasonal grazing to maintain adequate herbage and ground cover, (3) the 
stage of vegetation development to properly assess range readiness for 
grazing, (4) yearly effects of weather phenomena on forage development 
to plan cattle movements, (5) the effects of evapotranspiration and 
animal consumption on water ponds in grazing areas, and (6) the degree 
of range utilization in relation to its effects on long-term range 
productivity. 

This information if timely would help private range operators to 
assess animal-carrying capacities of large and inaccessible range 
holdings. Carneggie and DeGloria (1974) have shown that the repetitive 
coverage and synoptic view of LANDSAT can be a valuable tool for eval- 
uating rangeland resources in California. Deering et al. (1975) 
developed ,a green vegetation estimation model using LANDSAT data that can 
be used as an indicator of rangeland vegetation conditions in the Great 
Plains. Maxwell (1976) described how LANDSAT data can provide range 
managers with maps and tables giving standing herbaceous biomass for 
selected range types in northwestern Colorado. Thompson (1975) has 
shown that LANDSAT imagery can provide frequent estimates of rangelands 
in western Canada where estimate intervals now are 6 to 10 years apart. 
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Our objective was to investigate eight spectra? vegetation indices 
that can be determined from LANDS AT multispectral scanner (KSS) digital 
data (Richardson and Wiegand* 3977) to determine their usefulness to 
range managers for estimating animal harvestable forage in open rangeland 
areas of south Texas, 
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EXPERIMENTAL PROCEDURES 
Study Area and Field Methods 

The specific study area is located between 26°28 l and 26°42’ north 
latitude and 97°25’ and 97*49 ! west longitude, and includes about 81,000 
ha in Kenedy and Willacy Counties, in south Texas (Everitt et al., 1979). 
It is a transition zone between the Texas Coastal Prairies and the South 
Texas Plains vegetational regions (Gould, 1975). The Gulf of Mexico 
borders the area on the east. 

The topography is flat to gently sloping with elevations ranging 
from sea level to 15 m above sea level (U.S. Geological Survey 7.5 Minute 
Topographic Maps). 

The climate is mild, with short winters and relatively warm tempera- 
tures throughout the year. The average growing season exceeds 325 days 
(Texas Almanac, 1975). The 30 year average annual rainfall is 660 mm 
(Potter, 1976). Heaviest rains fall in May and September. 

Herbaceous production of primary interest to range managers relates 
to total standing herbaceous biomass available for grazing, without 
deleteriously affecting vegetation. Therefore, total standing herbaceous 
biomass samples were obtained, at or near the time of LANDSAT-2 clear 
overpasses fro® five different rangeland sites within the 81,000 ha study 

area. Since the LANDS AT MSS digital data is expected to be sensitive 

♦ 

only to the total standing green herbaceous biomass, it should yield best 
estimates of total standing herbage production when standing green herb- 
age predominates over standing brown herbage (Tucker, 1977). The five 
rangeland sites used were: 
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1. Tight sandy loam (improved grasses) . The brush had been controlled 
and the range reseeded with buffel grass ( Genchrus ciliaris L.). 
Native grasses have reestablished themselves. At the time of the 
study this was a highly productive grassland site. 

2. Tight sandy loam (allcia grass) . The brush had been controlled and 
the site has been improved by reseeding with alicia grass ( Cynodon 
spp , ) , At the time of the study, this was a highly productive grass- 
land site. 

3. Coastal sand (native grass ) . This native grassland site had few 
woody plants (less than 5%) other than an occasional mott of mesquite 
( Prosopis glandulosa Torr.) or live oak trees ( Quercus virginiana 
Mill.). At the time of the study, this was a moderately productive 
coastal grassland site. 

4. Deep sand (improved grasses) . The brush had been controlled and 
native grasses and herbs were reestablished. At the time of the 
study this was a moderately productive grassland site. 

5. Salty Flat (native grass ). A few salt tolerant herbaceous species, 
but no woody plant species were present at the time of the Study and 
productivity was low., 

We determined total standing herbaceous biomass production for each 
rangeland site by clipping all standing vegetation 2 cm above ground 
level in 20 quadrats (each 50 cm x 50 cm) (Stewart and Hutchins, 1936) 
once a month during the 1975 and 1976 range growing seasons. Percent 
canopy cover of woody vegetation was determined by the line-intercept 
method (Canfield, 1941). Besides these five rangeland sites in the study 
area, three other bare soil sites that were used are: 
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6. Tidal Flats . This site occurs on nearly level areas along the coast 
only a few feet above the mean high tide and has large bare soil 
areas or salt slicks. 

7. Cropland . These are blocks of cultivated land located in the 
southern part of the study area. Nonirrigated grain sorghum and 
cotton row crops are grown in the late spring and summer. 

8. Sand Dunes . These are coastal areas of deep unstable sand usually 
bare of vegetation, except for vegetation around their outer 
perimeter. 

These three bare soil sites (6 to 8) were not sampled for biomass 
measurements. For each of these sites, we assumed biomass was zero for 
regression analysis studies to relate biomass to LANDSAT spectral data. 
On July 31, 1976, the cropland site was planted to sorghum and so we did 
not use it for the regression analysis on this date. 

Although sandy mound live oak brush, deep sand mesquite bruoh , and 
tight sandy loam mixed brush rangeland sites were present in the study 
area (Everitt et al. , 1979), they were not used in regression analysis 
studies because their woody canopies obscured the herbaceous biomass . 
Lagunas were not used in the study because their signature varied from 
wet to dry seasons. 
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LANDSAT Digital Data 

We used MSS computer compatible tapes (CCT) and corresponding color 
images (Is 1,000, 000 scale) for four LANDSAT-2 overpasses. All four of 
the LANDSAT MSS bands were used, covering the 0.5- to 1.1-pro spectral 
region. These overpasses provided digital counts for a 185- by 185-km 
scene that included thG 81,000 ha rangeland study area in Kenedy and 
Willacy Counties in south Texas. Four cloud-free LANDSAT-2 overpass 
dates were chosen to relate satellite data to herbage production. An 
October 17, 1975 overpass (scene ID 2268-16190) provided an image of 
the area when most vegetation was in late season growth. A December 10, 
1975 overpass (scene ID 2322-16183) provided an image of the area when 
it was seasonally dormant \ near freezing air temperatures and a radia- 
tions! frost had occurred about a month earlier. The July 31, 1976 
overpass (scene ID 2556-16125) provided an image of the area when the 
vegetation was at peak growth after heavy rains in June and early July. 
The September 23, 1976 overpass (scene ID 2610-16112) provided a second 
data set comparable to the October 17, 1975 overpass. 

A supervised analysis was used with the MSS data that consisted of 
acquiring over 650 training pixels (picture element) from 15 training 
sites within the 81,000 ha rangeland study area (1.0% of total area) 
that intensive ground observations indicated were representative of the 
five rangeland and three bare soil sites previously described. The mean 
digital count for each spectral band of each site was used to calculate 
eight spectral vegetation index (VI) models, as described by Richardson 
and Wiegand (1977). Training sites were identified on gray scale maps 
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of the study area, from which record and pixel coordinates were deter- 
mined. The same training sites were used for all four LANDSAT overpass 
dates. Regression analysis was u£s<sd to test the ability of each of the 
eight VI models to estimate total herbaceous biomass production pre- 
viously determined for each of the five rangeland types. The regression 
model used is as follows: 

BIOMASS (kg/ha) = a Q * ^ VI, (1) 

where VI is any one of eight vegetation index models. 

Vegetation Index Models 

The same vegetation indices, used by Richardson and Wisgand (1977) 
to characterize the seasonal development of grain sorghum, were used in 
this study to quantify seasonal grassland changes. All models use 
LANDSAT digital data. A sun angle cosine correction factor (F) was 
used to correct for the seasonal sun zenith angle effects as follows 
(Deerin^ et al., 1975; Richardson and Wiegand, 1977): 

F = cos(? / cqs0, (2) 

where = 39°, tiedian sun zenith angle, and Q is the zenith angle of the 
LANDSAT overpass being processed. The raw LANDSAT CCT digital data were 
multiplied by F before the Vi's were calculated. 

The following formulas define the eight VI models: 

Transformed Vegetation Index Model (TVI), 

TVI = ( (MSS7-MSS5)/(MSS7+MSS5) t Q.5)\ (3) 

Transformed Vegetation Index Six Model (TVI6), 

TVI5 = ( ( MSS6-HSS5 ) /( MSS6+MSSS ) + 0.5) . (4) 
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Ratio Vegetation Index Model (RVI), 

RVI * MSS5 / MSS7. (5) 

Perpendicular Vegetation Index Model (PVI), 

PVI a (<Rgg5 - Rp5) 2 + (Rgg7 « Rp?) 2 )* 5 , (6) 

wher^i RggS a 0,851Rp5 + 0.35SRp7; Rgg7 a 0.355Rp5 + 0.148Rp7', RpS a MSS5; 
and Rp7 a MSS 7. 

Perpendicular Vegetation Index Six Model (PVI6), 

PVI6 a ((Rg g 5 - Rp5) 2 + (RggS - RpS) 2 )**. (7) 

where Rgg5 a -.438 + 0,543Rp5 t 0.498Rp6i Rgg6 = 2.734 + 0.498Rp5 J* 
0.457Rp6 ; Rp5 a MSSSj and Rp6 a MSSf? . 

Different^ Vegetation Index Model (DVI), 

DVI = 2.4MSS7 - MSS 5 . (8) 

Soil Brightness Index Modal (SBI), 

SBI = 0.433MSS4 + 0.632MSS5 + 0.586MSS6 + 0.264MSS7. (9) 

Green Vegetation Index Model (GVI), 

GVI = -0.290MSS4 - 0.562MSS5 + 0.600MSS6 t 0 ,491MSS7.(10) 

The TVI, TVI6 , and RVI models were used by Texas ASM University 
personnel as indicators of the amount and condition of rangeland vege- 
tation (Deering et al., 197S). The PVI, PVI6, and DVI models were 
developed by Richardson and Wiegand (1977) at Weslaco, Texas, for spec- 
tral monitoring of sorghum plant growth. The SBI and GVI models (Kauth 
and Thomas, 1976) were used in the Large Area Crop Inventory Experiment 
(LACIE) at the Johnson Spacecraft Center, Houston, Texas, for describing 
soil background reflectance and crop development, respectively. 
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RESULTS AND DISCUSSION 

Vegetation Indtx Models 

The results of the regression analysis that tested (equation 1) 
the ability of each of the eight VI models for estimating herbaceous 
biomas3 are given in Tables 1 through 4 along with the LANDSAT and 
ground observation data collected for each range site on each of the 
four sampling dates (October 17 and December 10, 1975; July 31 and 
September 23, 1976). All VI models yielded significant correlations 
with herbaceous biomass although the PVI and DVI models seemed to have 
generally higher correlation coefficients for all four dates. As ex- 
pected, the soil brightness index (SBI), a measure of soil background 
reflectance, was negatively (but insignificantly) correlated with biomass 
for any date. Since the individual LANDSAT band digital counts were not 
significantly correlated with biomass , the VI models are a considerable 
improvement for characterizing rangeland herbage production. 

Regression coefficients (r), relating total herbaceous biomass with 
PVI for each of the four LANDSAT overpass dates individually using open 
grass and bare soil range sites for all four dates, ranged from 0.74 to 

0. 97 (Table 5). The regression slopes and intercepts indicated that the 
biomass and PVI data follow about the same linear line for all four dates. 
Thus, the data for open range grass sites and bare soil for all four 
dates were combined into one regression, as shown in Table 5 and Figure 

1. The regression coefficient was significant (r = 0.90**) indicating 
that this single relation could be used to estimate standing herbaceous 
biomass production of open range from PVI at any time of the growing 
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season for either year even though the standard error of estimate, 722 
kg/ha, is larger than reported by Deering et alt (1977) and Harlan et. al* 
(1979) of 250 kg/ha for estimating natural vegetation. 

Phenological Variation of Rangeland Biomass 

The relative proportions of standing green, standing brown, and 
heads and stems biomass components constituting the total amount of 
biomass measured for the moderately productive coastal sand native grass 
rangeland site (site 3) for the 1976 growing season are shown in Figure 
2, Standing green biomass accounted for the large increase in total 
biomass starting in June and ending about October. The increase in 
heads and stem biomass from March through October, with an inexplainable 
dip in September, is consistent with the known phenology of grass growth. 

amount of biomass for standing brown was unexplainably variable 
throughout the growing season. 

The seasonal changes in amounts of total herbaceous bionass for the 
moderately productive coastal sand native grass (site 3) and the less 
productive deep sand improved grass (site 4) rangeland sites are shown in 
Figure 3. Total herbaceous biomass for both sites was about the same until 
July and August, when the coastal sand site produced more green biomass, 
and therefore, more total biomass than did the deep sand site. 
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SUMMARY AND CONCLUSIONS 

Our results showed that LANDSAT data are significantly correlated 
with total standing herbaceous biomass. In the study it was not possible 
to test the combined ability of LANDSAT and weather data to estimate 
biomass conditions, because cloudy weather prevented LANDSAT coverage 
on all, except 2 of the 8 biomass collection dates, in both 1975 and 
1976 resulting in insufficient data. However., weather information has 
been shown to be a good estimator of herbaceous biomass (Deering et al., 
1975), thus it seems likely that the combined use of LANDSAT and weather 
data would improve the reliability of herbaceous biomass production 
estimates . 

Everitt et al. (1979) showed that it is possible to distinguish 
among open grasslands, mixed brush, and live oak rangeland types. Thus, 
once the areal extent of open grasslands have been remotely inventoried, 
using LANDSAT MSS data-; then the relations developed in this paper could 
be used to estimate the herbage production of these areas. For areas 
that are largely open grassland, like those studied by Carneggia and 
DeGloria (1974) and Maxwell (1976), these procedures could be used 
directly. Therefore, it seems likely that this information could be 
displayed as computer generated maps and tables for use by range managers 
for estimating herbage production of open rangeland. 
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Table 5. Regression analysis of total standing herbaceous biomass 
(BIOMASS) with the perpendicular vegetation index (PVI) 
determined from LANDSAT MSS bands 5 and 7 for four LANDS AT 
overpass dates individually and all dates combined. 


LANDSAT 




OVERPASS 

LINEAR REGRESSION EQUATION 

r 

Sy.x 

DATE 



(kg/ha) 

10/17/75 

BIOMASS = -370 + 277 PVI 

0,97 

527 

12/10/75 

BIOMASS = -71 + 186 PVI 

0.74 

684 

7/31/76 

BIOMASS = -180 + 250 PVI 

0.87 

1063 

9/23/76 

BIOMASS = -264 + 231 PVI 

0.86 

812 

ALL DATES 





BIOMASS = -297 + 253 PVI 

0.90 

722 

COMBINED 
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FIGURE CAPTIONS 

Figure 1. Overall regression (solid line) of total standing herbaceous 
biomass collected on October 17 and December 10, 1975, and 
July 31 and September 23, 1976, with the perpendicular 
vegetation index determined from LANDSAT MSS bands 5 and 7. 

Figure 2. Temporal plot of total, standing green, standing brown, and 
heads and stem herbaceous biomass collected for a coastal 
sand native grass range site , located in south Texas , for 
the 1976 growing season. 

Figure 3. Temporal plot of total standing herbaceous biomass for medium 
(coastal sand native grass) and low (deep sand improved grass) 
productive forage range sites in south Texas for the 1976 
growing season. 
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